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magnitude but opposite sign. To accomplish this, we compute the relevant components of the 
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the full quark propagator, obtained as a solution of the quark gap equation, and the gluon propa- 
gator and ghost dressing function, obtained from large-volume lattice simulations. 
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1. Introduction 

In the last few years, considerable progress has been made in our understanding of the infrared 
(IR) behavior of the fundamental Green's functions of QCD, such as gluon, ghost, and quark prop- 
agators ©g^B^BSBil' 88 wdl as some of the basic vertices of the theory |g 0, 011], 
and their relation to the confinement and dynamical chiral symmetry breaking (CSB) [|[]. In fact, 
there is a broad consensus that one of the most important ingredients for the CSB is the non-abelian 
quark-gluon vertex, which controls the way the ghost sector enters into the gap equation. Specif- 
ically, this vertex introduces a numerically crucial dependence on the ghost dressing function and 



the quark-ghost scattering amplitude [|13|]. This latter quantity satisfies its own dynamical equation, 
which may be decomposed into individual integral equations for its various form factors. Here we 
will present the first steps towards the determination of the longitudinal quark-gluon vertex form 
factors for a particular kinematic configuration: the quark symmetric limit where the incoming and 
outgoing quark momenta have the same magnitude and opposite signs. To do that, we compute 
numerically the relevant quark-ghost scattering kernel components at the "one-loop dressed" ap- 
proximation, at the same kinematic point, using as ingredients the nonperturbative lattice results 
for the gluon propagator and ghost dressing function of Ref. [0], and the solution of the quark gap 



equation obtained in Q13| ] for the full quark propagator. 



2. Ingredients and definitions 

Consider the conventional quark gluon vertex shown in Fig. [I], and defined according to 
i^ Vj xif,(q,P2,-pi) =igtf j r tl (q,p 2 ,-pi); r$\q,p 2 ,-pi) = Y tl ; q+Pz=Pi- (2.1) 




Figure 1: The conventional quark-gluon vertex with the momenta routing used throughout the text. 



In the Batalin-Vilkoviski (BV) formalism, the Slavnov-Taylor identity (STI) satisfied by this 
vertex reads (in the kinematic configuration chosen) [|]] 

q^A% Vj Wi(q^P2, ~Pi) = F{q 2 ) Ty^p^Ty.^^q, -p\) +T^ w . ca (p 2 , -puq)Ty,y k {pi) 

(2.2) 

In the formula above F(q 2 ) denotes the so-called ghost dressing function which is related to the 
full ghost propagator D ab (q 2 ) through 



iD ab (q 2 ) = i8 ab ^-; F W(q 2 ) = \, 



(2.3) 
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while is the inverse of the full fermion propagator Sjj(p) obtained by solving the equation 

iS u (p)T wm = 8 ik ; isf (p) = iSijS® (p) = iS u . (2.4) 

The standard decomposition for the inverse of the full quark propagator S^ 1 (p) is given by 

iS- l (p)=i[A(p 2 )0-B(p 2 )], (2.5) 

where A(p 2 ) and B(p 2 ) are, respectively, the Dirac vector and scalar components. 

In addition, and represent the so-called antifields associated to the spinor fields \j7 k and 
Yk respectively; they have ghost charge -1, (mass) dimension 5/2, and obey Bose statistics. The 
Green's functions r v *y jC a and T^. d r^> are shown in Fig. ||. 

These two functions are not independent, being related by "conjugation"; indeed, to get one 
from the other, we need to perform the following operations: (i) exchange —p\ with p2- —p\ P2', 
(z'z) reverse the sign of all external momenta: q,—p\,p2 —q,p\,—p2\ {Hi) take the hermitian 
conjugate of the resulting amplitude. 

Then, introducing the function 

H ij{q,P2,-p\) =gtfjH(q,p 2 ,-pi) = -iT VjC ay*(p 2 ,q,-pi), (2.6) 
H"j(-q,Pi,-p2) =gti' j H(-q,p l ,-p 2 ) =ir v * WiCa (p 2 ,-pi,q), 



and factoring out a the common color and gauge coupling combination gtf=, the STI ( |2.2| ) can be 
rewritten as 

q^r^q,p2,-pi)=F(q 2 )[S- l (p l )H(q,p 2 ,-pi)-H(-q,Pu-P2)S-\p2)], (2.7) 

with H obtained from H through the set of operations detailed above. The H function admits the 
general form factor decomposition [ p"4| 

H{q,p 2 ,-pi) =XbI+Xi>i+X 2 > 2 +X 3 cf /lv ^, (2.8) 

where the form factors X t are functions of the momenta, X,- = Xj(q 2 , p 2 ., p 2 ) and o^ v = 1 /2[y^ , y v ] 
(notice the i difference with respect to the conventional definition of this quantity). One then obtains 
automatically the expansion 

H(-q,p h -p 2 ) =X I+X 2 M^MWiivPiP2, (2-9) 

where X =Xi(q 2 ,p\,p\). 

At tree-level, one clearly has X (0) = X^ = 1, with the remaining form factors vanishing. 

The most general Lorentz decomposition for the longitudinal part of the vertex F^(q,p2, —p\) 
appealing in Eq. (2.7) can be written as fll-i] ] 

Fn(q,P2,-pi) =L\Yn +h)(Pi+P2)n +L3(Pl+P2)n+ L 4°Liv(Pl+P2) V , (2.10) 

where L, are the form factors, whose dependence on the momenta has been suppressed, in order to 
keep a compact notation, i.e., L, = Li(q 2 ,p 2 ,p 2 ). Notice that the tree level expression for Tn is 
recovered setting L\ = 1 and L2 = L3 = L4 = 0; then, = y^. 
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Figure 2: The auxiliary functions through which the STI satisfied by the quark-gluon vertex is satisfied. For 
convenience we show the momenta routing matching the kinematics chosen for the quark-gluon vertex as 
well as the tree-level and one-loop terms in the corresponding loop expansion. 



Contracting Eq. ( |2.10| ) with , we have 

= {pi - pl)L 3 I + [{p\ - pi)L 2 + L\]$\ + [(p\ - p\)L 2 - Li]£ + 2L 4 a^pl (2. 1 1) 
In addition, substituting into Eq. ( [2^ ) the full quark propagator S~ l (p) of Eq. (2^5), and the 



expressions for H and H given by Eqs. ( |2.8D and ( |2.9| ) respectively, we find that the rhs of Eq. 1 1| ) 
can be also expressed in terms of the functions A, B and Xj's. Then, it is relatively straightforward 



to demonstrate that the L,'s may be expressed as [13] 



+ 



2 

F{q) 



{A(p l )[X -(pj + p v p 2 )X 3 }+A(p 2 )\X -(p 2 2 + pvP2)X 3 }} 
{B( Pl )(X 2 -Xj) + B(p 2 )(X 2 -Xj)}; 



L 2 = *, F 2 9) 2, {A(p i )[Xo + (pj-p l -P2)X 3 ]-A(p 2 )\X + (p 2 2 - PvP2 )X 3 }} 



{5(/>i)(Xi +X 2 )-B(p 2 )(X l + X 2 )} ; 



L 3 = {A(pi) (p\x x +p vP2 X 2 ) -A{ P2 ) {p 2 2 Xi+pvp 2 X 2 ) -B( Pl )X + B{p 2 )X } ; 

P1-P2 



L 4 



{A( Pl )X 2 -A(p 2 )X 2 -B{p x )X 3 +B(p 2 )X 3 } . 



(2.12) 



It is interesting to notice that setting in Eq. ( |2.12| ) Xq = Xo = 1 and X,- = X; = 0, for i > 1, and 
= 1, we obtain the following expressions 



A(pi)+A(p 2 ) A( Pl )-A(p 2 ) B{p x )-B( P2 ) 
L \ — ~ > L 2 — n , 7 t\ — 1 



2(P?-P 2 2 ) ' 



2 2 
P\~P 2 



L 4 = 0. (2.13) 



which give rise to the so-called Ball-Chiu (BC) vertex [15], widely employed in the literature for 
studies ofCSB [& 
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3. The "one-loop dressed" approximation for H 



It is clear from Eq. (2.12), that in order to determine the longitudinal form factors L, , it is 
necessary to know the nonperturbative behavior of the form factors X,-. 

To obtain a nonperturbative estimate for H and its form factors, we will study the "one-loop 
dressed" contribution represented in the diagram of Fig. |[ and given by 



# [ %P2,-Pl) = 1 " \ 




Figure 3: The quark-ghost scattering kernel at "one-loop dressed" approximation. 



#W (q,p 2 ,-pi ) = 1 - \ iC A g 2 J^ v (k)G v (k - q)D{k - q)S{k + p 2 )T il {k,p 2 ,-k - p 2 ) , (3. 1) 

where Ca is the eigenvalue of the Casimir operator in the adjoint representation, A^ v (q) is the full 
gluon propagator; in the Landau gauge 



<v(?) 



id ab 



A(q). 



(3.2) 



For evaluating Eq. ([O]) further, we will use the following approximations: (i) the full gluon- 
ghost vertex will be replaced by its tree-level value, i.e. iGy hc = —gf abc (k — q) v . Note that, since 
the full A^ v (k) is transverse, only the q v part of the gluon-ghost survives since k v A^ v (k) = 0; (ii) 
for the vertex we will use the following Ansatz 



rJk,p2,-k-p 2 ) 



+ 



2 

F(k) W 



[A{p 2 + k)+A{p 2 )]y^ 

[A{p 2 + *) -A(p 2 )} (l^ + k 1 )-! [B(p 2 + k)- B( P2 )} 



.(3.3) 



Notice that the above Ansatz satisfies the STI of Eq. ( |2.7D when H = 1. Again, due to the 
transversality of A^ v (k), the second term on the rhs of Eq. (3.3), which is proportional to the 



longitudinal momentum k^, does not contribute in the Eq. ( |3.1| ). 



Then, inserting into Eq. ( |3.1[ ) the propagators of Eq. ( |3.2[ ) and the Ansatz for the quark-gluon 

ve the following expre 

,Pz,~Pi) = 1 + /2 -n 1 / ^iP2,q,k)fip2,q,k), (3.4) 



vertex given by Eq. ( |3.3| ), it is straightforward to derive the following expression for H 

H [i],. _ _ n i , .8 2 Ca 
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where 



JtT(p 2 ,q,k) 



D{k - q)F(k)A(k) [A(p 2 + k) +A(p 2 )} 



A 2 (p 2 + k)(p2 + k) 2 -B 2 (p 2 + k) ' 
while all spinorial structure is included in 



f(j>2,q,k) =A(p 2 + k) 



q ■ k 

h4+$4-k-q-fc$-jp 



+ B(p 2 + k) 



k^- 

* k 2 



(3.5) 



(3.6) 



4. Quark symmetric configuration 



The projection of the form factors X, appearing in the definition (2.8), for arbitrary kine- 
matics boils down to a complicated system of several equations. In order to make the problem 
at hand technically more tractable, we will only compute it in a specific kinematical limit: the 
quark symmetric limit where the quark momenta have the same magnitude and opposite signs, i.e. 
Pi = -p/2, p 2 = p/2 and q = -p . 

In this kinematical configuration, it is easy to see that only the (Dirac) form factor L\ survives, 
and the vertex of Eq. ( 2.10| ) simplifies to 



r tl (-p,p/2,p/2)=L 1 y fl ; where L\ = L\(p 2 ,p 2 /A,p 2 /A). 



(4.1) 



In addition, the quark-ghost scattering kernel H and H, given by Eq. ( |2.8| ) and Eq. ( |2.9[ ), simplify 
in this limit, and the terms proportional to p 1 become linearly dependent such that 



H(-p,p/2,p/2) = X Q l+(X 2 -X l )p / /2, 
H{p-p/2-p/2) = X I + (X l -X 2 )^/2. 



(4.2) 



Setting p\ = p 2 = p/2 and q = —p in Eq. (p.4[), and taking the appropriate traces, it is straight- 
forward to derive the following expressions for the form factor Xo and the subtraction X2 — Xi (in 
the Euclidean space) 



X (p) = 1 + 



g 2 C A rD{p + k)F{k)\{k)A 2 [A 2 +A { 



A 2 2 (p/2 + k) 2 + Bj 



(p-k) 



21 



k 2 



X 2 (p)-X l (p) 



g 2 C A fD(p + k)F(k)A(k)B 2 [A 2 +A l 
2p 2 



A 2 {p/2 + k) 2 +B 2 



{p-kf 



k 2 



(4.3) 



where A 1 = A(p/2), A 2 = A(p/2 + k), Bi =B(p/2), and B 2 = B(p/2 + k). Due to the fact that the 
momenta p\ and p 2 have the same magnitude in this configuration, it is possible to show from the 
definition (ph that X = X and X 2 -X l =X 2 -X l . 



Then, it is straightforward to see, from Eq. ( p.l2[ ), that in the symmetric quark limit L\ simpli- 
fies to 



F(p)X Q (p)A(p/2)+F(p)[X 2 (p) -Xi(p)]fl(p/2) . 



(4.4) 



As we have seen in Eq. (4.3), both form factors Xo and X2 — X\ depend on the nonperturbative 
form of the four basic Green's functions, namely A(q), F(q), A{p) and B(p). Therefore, in order 
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Quark-ghost scattering kernel 

Scalar component 

Quark symmetric configuration 



1E-3 0.01 




Quark -ghost scattering kernel 
Dirac components 
Quark symmetric configuration 



P [Gen 



p [GeV ] 



Figure 4: The scalar form factor Xq (left panel) and the combination of the form factors X2 — X\ (right panel) 
in the symmetric quark configuration for <x(jJ, 2 ) = g 2 /4-n = 0.30. 



to proceed with the numerical analysis, we use for A(g) and F(q) the lattice data obtained by 
[Q], while for A(p) and B{p) we use the solution of the quark gap equation obtained in Ref 
All these functions were renormalized at pL = 4.3 GeV, and in all our calculations we have fixed 
a^ 2 ) = g 2 /47r = 0.30. 

In Fig. I], we show the numerical results for Xo(p) (left panel), and X2(p) —X\ (p) (right panel). 
On the left panel of Fig. |], we can see that Xo shows a maximum located in the intermediate 
momentum region (around 1 — 2 GeV ), while in the UV and IR regions the curve goes to its 
perturbative value i.e. Xo — > 1. 

On the right panel of Fig. ^| we notice that the combination X2 — X\ saturates at a finite value 
in the deep IR region, while in the UV it vanishes asymptotically. 

With all ingredients available, we are now in position to determine the behavior of the Dirac 
form factor L\ in the symmetric quark configuration. 

In Fig. H we show the result for L\ obtained from Eq. (4^4). As we can clearly see, LI develops 
a sizable plateau in the IR region, and as expected, it recovers its perturbative value in the deep UV 
region. 



Form factor - Vertex 

Quark symmetric configuration 



1E-3 0.01 0.1 1 10 100 

p 2 [GeV 2 ] 



Figure 5: Numerical result for the vertex form factor L\ in the quark symmetric configuration when 
a(ju 2 ) = 0.3. 
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5. Conclusions 

We have presented the general methodology for determining the longitudinal form factors of 
the quark-gluon vertex from the fundamental STI that it satisfies. A key ingredient in this analysis 
is the quark-ghost scattering kernel, H and its conjugate H, whose field-theoretic origin and basic 
kinematic properties are rather subtle. 

The first nonperturbative estimate of the form factors comprising H has been computed using 
the "one-loop dressed" approximation of the corresponding integral equation, under certain rea- 
sonable dynamical assumptions. For the purposes of this presentation we have limited our analysis 
to the particular kinematic limit known as "quark symmetric configuration", which gives rise to 
considerable technical simplifications. The Dirac form factor of the quark-gluon vertex, L\, has 
been obtained in this particular kinematic configuration, and its basic features have been studied. 
The methodology presented here may be directly extended to arbitrary kinematic configurations, 
furnishing valuable information on such a fundamental quantity as the quark-gluon vertex. 
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